
set scheme lean2
use "$altered_data/event_study_vaccination.dta", clear

//Days relative to end of summer break
collapse cases_melde, by (age rel_time)
xtset age rel_time


twoway (tsline cases_melde if age ==1)(tsline cases_melde if age ==2)(tsline cases_melde if age ==3)(tsline cases_melde if age ==4), ///
		title("`header'", color(black)) ylabel(0(15)60, nogrid angle(horizontal)) ///
		graphregion(color(white)) ytitle("Daily cases per 100k") xtitle("Days relative to end of summer break") ///
		legend(order(1 2 3 4) label(1 "Age: 5 - 14") label(2 "Age: 15 - 34") label(3 "Age: 35 - 59") label(4 "Age: 60+") cols(2) pos(6)) xlab(-15(15)45) 
graph export "${plots}incidence_openings.png", replace width(3000)
	
		
//Days relative to start of summer break
use "$altered_data/event_study_vaccination.dta", clear

collapse cases_melde, by (age rel_time_c)
xtset age rel_time


twoway (tsline cases_melde if age ==1)(tsline cases_melde if age ==2)(tsline cases_melde if age ==3)(tsline cases_melde if age ==4), ///
		title("`header'", color(black)) ylabel(0(15)60, nogrid angle(horizontal)) ///
		graphregion(color(white)) ytitle("Daily cases per 100k") xtitle("Days relative to start of summer break") ///
		legend(order(1 2 3 4) label(1 "Age: 5 - 14") label(2 "Age: 15 - 34") label(3 "Age: 35 - 59") label(4 "Age: 60+") cols(2) pos(6)) xlab(-15(15)45) 
graph export "${plots}incidence_closures.png", replace width(3000)
